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The possibihty of facing critical phenomena in nuclear fragmentation is a topic of great interest. 
Different observables have been proposed to identify such a behavior in particular, some 

related to the use of information entropy as a possible signal of critical behavior M. In this work 
we critically examine some of the most widespread used ones comparing its performance in bond 
percolation and in the analysis of fragmenting Lennard Jones Drops. 
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When nuclei collide at energies of the order of hundreds of MeVs, the resulting highly excited system breaks up 
into several medium-size fragments. This phenomenon is usually referred to as multifragmentation and attracts a 
lot of attention in the nuclear community because it is conjectured that the fragmenting system might undergo a 
phase transition. In particular the possibility of facing critical behavior (originally triggered by the seminal work of 
] the Purdue Group [^|) has induced the development of different signatures to characterize such a behavior from the 
1^ . analysis of the observed mass spectra. Starting with the work of Bauer |^ and Campi [|j the percolation model has 
OO [ emerged as a quite useful test-bed for the analysis of the tools developed to study nuclear multifragmentation. 
04 , In 1^ the information entropy (Si) has been proposed as a signature of critical behavior and two systems has been 
analyzed: a lattice gas model and a molecular dynamics model. 

In this work we test the information entropy together with other proposed signatures of critical behavior in a system 
(percolation on a simple cubic lattice) that it is known to undergo a second order phase transition. For this model, 
^ the behavior near the critical point (including the critical exponents) is well known. Besides the information entropy, 
we also analyze the normalized variance of the size of the maximum fragment (NVM) the second moment of the 
T— I distribution (M2) the best fitting power law ||^ and minimum t 

We have calculated the above mentioned signatures from the mass distributions resulting from a bond percolation 

■ process on two simple cubic lattices of sizes 30 x 30 x 30 and 6x6x6. It is known that for the larger one, the finite 
, size effects are negligible, meanwhile for the smaller one (whose size is of the order of the one corresponding to nuclear 

■ system's size) the finite size effects are important 

"V* '. Once the various criticality signatures have been tested we use them to analyze excited Lennard-Jones drops 
' Q ■ undergoing fragmentation. We focus on whether the non-equilibrium process of fragmentation can be characterized 
as a critical phenomenon. 

^ ■ In percolation, the number of clusters of size A per lattice site, ua, at a given bond activation probability, p, was 
proposed in M to be given by : 



UAip) = qoA-^fiz), (1) 

where qo is a normalization factor, r a critical exponent and f{z) a scaling function. This equation shows that 
nA{p) / qoA~'^ depends onp and A via the combination z = {{p — pc)/pc)A'^ = eA'^. Thus it is clear that at the critical 
point n^(pc) = 9o^~^/(0) with /(O) = 1. All this is valid for not too small fragments and infinite systems for which 
the percolating cluster has been removed. In finite systems we won't get a pure power law because finite size effects 
introduce distortions. Nevertheless, there exists a range of masses (depending of the size of the system) for which the 
mass spectra presents a power law like behavior, as can be seen in figure (If). 

There are two main methods to find the critical point from the mass spectra. The first one consists in fitting the 
mass spectra with a simple power law, with two independent parameters, qq and t, and then look for the value of p 
(or the quantity taken as control parameter) for which the exponent r is minimum |^,^ . 

The second one was developed in |^ and takes care of the constraint imposed by the normalization, i.e. qo ~ qa(T). 
This dependence can be found by analyzing the first moment of the cluster distribution at the critical point |]1G|| : 

Mi(e = 0) = ^ n^(e = 0)A = go ^ A'-^ = 1, (2) 

A A 

with the sum running over all clusters. Equation (^ shows that the normalization constant, qo, depends on r via 
a Riemann C function: 
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Therefore, if we take into account the dependence of go with r via equation (H), we have to fit the cluster distribution 
with a single parameter power law. In this way, the critical point corresponds to the best fitted spectra, i.e. minimum 

The critical point can also be found by looking for the maximum (divergence in an infinite system) of the moments 
of the distribution [Q or of the normalized variance of the size of the maximum fragment (NVM). 
The moments of the cluster distribution, Mk, are defined as: 

Mfc=^A'=nA, (4) 

A 

where ua is the cluster distribution. 

The normalized variance of the size of the maximum fragment (NVM) is defined as: 

NVM < ^ma.- < Ama. 

^ ^max ^ 

where Amax is the largest fragment of a given event and < ... > is an average over an ensemble of events at a given 
probability p. 

Finally, the information entropy, as proposed in is defined as: 

Si = -^PaMpa), (6) 
A 

where pA is the probability of detecting a fragment of mass A defined as pA = NA/Nt, with Na the number of 
fragments of size A detected in a set of experiments and Nt — J2a the total number of fragments detected in the 
same set of experiments. This magnitude is referred to as 5*1 because it belongs to the family of generalized Renyi 
entropies defined by pl| ]: 

In figure (1) we show the results of our calculations for percolation three dimensional lattices of two different sizes 
specified above. The NVM in figure (la) , M2 in figure (lb), in figure (Ic) (the minimum of denotes the best 
fitting power law with one-parameter fit), the r^ff in figure (Id) (corresponding to a two-parameter fit of the cluster 
distribution. The minimum value of r would correspond to the critical point) and S'l in figure (le). 

We can see that for the mass spectra corresponding to a lattice of size 30 x 30 x 30 (filled symbols of figure (1)) 
the signals NVM, M2 and give the same critical probability {pc = 0.25 ± 0.01). The extremes are very sharp 
and the critical probability is the one corresponding to the infinite system limit, indicating that the finite size effects 
are negligible. On the other hand, the maximum of the information entropy is not as sharp and it is also shifted (at 
p = 0.225 ± 0.015) with respect to the extreme of the other signatures, showing that it is not a good signature of the 
critical point. Finally, we can see that when we fit the cluster distribution with a power law with two independent 
parameters the exponent r is minimum at p = 0.24 ±0.01, slightly shifted from the other signals (but with overlapping 
error bars) and with a minimum less pronounced than when we fit the mass spectra with r and qo{T) (minimum x^)- 

For the smaller (6x6x6) lattices (open symbols in figure (1)), we see that the peaks of the various signatures 
analyzed are slightly shifted from one another but with overlapping error bars, even for the case of the information 
entropy. The best fitted mass spectra according to the condition given by equation (||) is at p = 0.31 ± 0.01 (figure 
(Ic), the NVM is maximum at p = 0.30 ± 0.01 (figure (la)), the M2 at p = 0.29 ± 0.01 (figure (lb)) and Sl{p^) at 
p., = 0.27 ±0.05 (figure(le)). 

As we can see, when the system is small, finite size effects can introduce distortion and smear out the signature of 
the critical point. On the other hand, 5*1 is shifted with respect the critical point when finite size effects are negligible. 

In order to understand the behavior of the information entropy in this kind of systems, we have chosen to generate 
a family of clusters distributions according to equation (|l|) with the critical exponents corresponding to percolation 
in an infinite lattice in three dimensions. In this case, the normalization constant qo is taken as the inverse of the 
zeta Riemann function evaluated at r = 2.18 and the scaling function f{z) is the one fitted in |l^. We have also 
chosen pc = 0.25 i.e. the infinite size value. These cluster distributions have the same functional form as the ones 
corresponding to a percolation problem, but without the percolating cluster and without finite size effects. 
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The information entropy can then be calculated by replacing nA{p) in the definition of Si (equation @) by its 
expression according to equation (jl|) . The results of this calculation are plotted in figure (2c) and we can see that 
the maximum in the entropy does not correspond to the pure power law distribution (pc = 0.25 in this case) as was 
claimed in [|| but to the distribution corresponding to p = 0.22 showed in figure (2d). We can check this result by 
calculating the derivative of against e = {p — Pc)/pc- 

^ - 90 E 4"'^4"'^^?/'(^i, ^)/(^2, e)log[{s^-fisi,r))/isrf{s,,r))], (8) 

SlS2 

with /(s, e) — f{z) the scaling function and f'{z) = If we evaluate equation (^ at the critical point, i.e. e = 0, 
we get: 

= 727;) ^'C^) E s"is[-^\si-^\logisi) - logis2)] ^ -1.6. (9) 

^ ' S1S2 

It is clear from this expression that the information entropy does not show an extreme in the critical point for this 
family of distribution functions and, therefore, 5*1 is not a signal of critical behavior. 

In Figure (2a) we also plot the coefficient for the cluster distribution fitted with one parameter (go — Qo{t) 
according to equation (^)) and in figure (2b), the second moment of the distribution, M2. It can be seen that these 
magnitudes signal the critical point very accurately. 

Finally, we apply these tools to molecular dynamics simulations of fragmentation. These simulations have been 
described in previous works (see p3|-p^ for details). We follow the time evolution of excited drops composed of 147 
particles interacting via a Lennard-Jones potential that undergo multifragmentation. 

Using this model and taking the excitation energy as a control parameter we perform the same analysis as in the 
case of the percolation model. The corresponding results are displayed in figure (3). In this figure we can see that 
M2 and NVM peak at £' = 0.3 ± 0.2e, meanwhile the best fitted mass spectra according to equation (||) (minimum 
X^) is, approximately, ai E — 0.15 ± 0.15e. On the other hand the information entropy monotonically increases in 
the range of energies analyzed and therefore does not show a maximum when the system displays a power law mass 
spectra, i.e. in the energies around E = 0.2e. 

In conclusion, we have analyzed different quantities proposed in the literature as signatures of critical behavior: 
i) the information entropy (Si), ii) the normalized variance of the size of the maximum fragment (NVM) , iii) the 
second moment of the distribution (M2), iv) the best fitting power law and minimum t . We found that 5*1 is not 
a reliable signature of a second order phase transition as was proposed in For this purpose, we have analyzed 
systems that display critical behavior: percolation on lattices of two sizes. For the larger one, where the finite size 
effects can be neglected, it is clear that the signature corresponding to information entropy is shifted with respect 
to the other quantities analyzed. But when the size is much smaller and the finite size can not be neglected, the 
signatures analyzed in this work display peaks that do not coincide in a given value of the probability. Instead, they 
are smeared in a small range (see figures (1)). In this case, the broad peak of the entropy can turn out to be rather 
close to this range of values (even though the error is large) and then the information entropy can be confused with 
the signatures. 

This behavior is confirmed in the analysis of highly excited Lennard-Jones drops that undergo multifragmentation. 

Finally, according to these results, we would choose as signatures of critical behavior the following signals: the best 
fit mass spectra with a one-parameter power law (minimum x^), the normalized variance of the size of the maximum 
fragment (NVAl) or the second moment of the cluster distribution (M2) and we would discard the information entropy 
iSi). 
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Figure Captions 

Fig.l: The normalize variance of the size of the maximum fragment, NVM (a), the second moment of the cluster 
distribution M2 (b), the coefficient of a one parameter power law fitting (with qo — qoir)) (c), the r exponent 
from the fitting of two parameters power law (d) and the information entropy (6*1) (e) against bond probability for 
percolation lattices of sizes 30 x 30 x 30 (full triangles) and 6x6x6 (empty circles). In (f) we plot the cluster 
distribution for the critical probabilities at both sizes. 

Fig. 2: The coefficient (a) and M2 (b) against probability for the cluster distributions generated from equation 
(|l|) using the scaling function f{z) obtained in fl^ . In (c) we plot the information entropy calculated by replacing 
the cluster distribution given by eq.(|l|) in the definition of Si (equation (^)). In (d) we can see both: the cluster 
distribution corresponding to the critical probability (chosen at pi = 0.25, a pure power law) and the one that 
maximize the information entropy {pi — 0.22). Here it is clear that the distribution that maximize the entropy is not 
a power law. 

Fig.3: The coefficient (a), the M2 (b), the NVM (c) and 5*1 against the energy per particle for a molecular 
dynamics multifragmentation process for 147 particles interacting via a Lennard- Jones potential. 
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